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I Abstract 

I I Initial data for boosted Kerr black hole are constructed in an axi- 

^ . ally symmetric case. Momentum and hamiltonian constraints are solved 

ijr^ ' numerically using finite element method (FEM) algorithms. Both Bowen- 

fT^ , York and puncture boundary conditions are adopted and appropriate re- 

T-H . suits are compared. Past and future apparent horizons are also found 

' numerically and the Penrose inequality is tested in detail. 

o : 

^ PACS numbers: 04.20. Cr, 04.20.Dw, 04.70.Bw 

o : 

qhI 1 Introduction 

gj-jl We will construct initial data for single boosted Kerr black hole in the axially 

• • . symmetric case with the total momentum parallel to the angular momentum of 

. ^ ' the black hole. We assume also asymptotical flatness. Two popular approaches 

. to the constructions of black hole initial data will be adopted. The first is 

5h ' the Bowen-York solution jlj based on the conformal transverse-traceless decom- 

position [3] of Einstein equations with certain additional assumptions such as 
conformal flatness and maximal slicing. The second is the puncture approach 
12] developed by Brandt and Briigmann. Their idea was to remove analytically 
the coordinate singularity at the location of the black hole. In this way they 
avoided inner boundary conditions present in the Bowen-York approach. We 
follow these ideas as close as possible but some important modifications are 
necessary. Garat and Price [H] have proved nonexistence of conformally flat 
slices of the Kerr spacetime. Thus the assumption of conformal flatness must 
be abandoned. We assume maximal slicing in order to simplify the problem 
but in general it is not necessary. Our approach is straightforward. We use 
the quasi- isotropic radial coordinate |3] , 0] in which it is easy to show that the 
Kerr solution is composed of two isometric regions smoothly joined at a sphere 
of some radius. Then the original Bowen-York inner boundary conditions can be 
applied. The momentum constraints in the axially symmetric case can be solved 
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analytically for the flat background metric [^S], [201 j and they can be solved nu- 
merically in the case of metrics which are not conformally flat. We solve the 
hamiltonian constraint using the finite element method (FEM) [5] techniques, 
then we find the apparent horizon and test the Penrose inequality. 

Our paper consists of three parts. In the first we shortly review Einstein 
constraint equations. The second section discusses the numerical approach and 
special assumptions that are made there. The third part deals with detailed 
results. 

2 Einstein Constraint Equations 
2.1 ADM framework 

In the Cauchy formulation of Einstein equations [HliISI the whole 4-dimensional 
manifold is foliated into a set of the so-called slices. These slices are the 3- 
dimensional, spacelike surfaces labeled by a parameter t (time). The Einstein 
equations are projected onto these hypersurfaces. One obtains six evolution 
equations and four constraint equations that the metric, 7^, and the extrinsic 
curvature, Kij, of each slice must satisfy. In the ADM (3-fl) decomposition of 
the spacetime the metric can be written as 

ds^ = -N^dt^ + 7y (da;* + l3'dt){dx^ + (3^dt). (1) 

The proper time between slices is given by the lapse function N and the co- 
ordinate drift between slices is described by the shift vector /3\ The evolution 
equations in a vacuum have the following form 

idt-Lp)j,j = -2NK,j, (2) 
{dt-Lp)K,, ^ {-W,Wj+R,,+KK,,~2K,kK!^)N. (3) 

Here Rij is the Ricci tensor on the slice, K is the trace of the extrinsic curvature 
and the covariant derivative Vi is taken with respect to the spatial metric jij. 
Lp denotes the Lie derivative along the shift vector /3'. 

The minimal set of initial data consisting of the 3-dimensional metric 7^ 
and the extrinsic curvature Kij cannot be specified independently on the initial 
slice. These data are constrained by four equations, as mentioned above, because 
the slices must fit properly into the 4-dimensional manifold. The first of these 
equations is known as the hamiltonian (or scalar) constraint 

R + - K,jK'^ ^0, (4) 

and the other three, called momentum (or vector) constraints read as 

Wj(K'^ --f'^Kl) = 0. (5) 

Here R denotes the Ricci scalar. Let us note that if the above constraint equa- 
tions are satisfied on one initial hypersurface, the same holds true on each slice. 
It is also important that the constraint equations do not depend on the lapse 
function N and the shift vector /3* which can be freely chosen on each slice. 
This is the base of the (general) gauge invariance of the theory. 
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2.2 York-Lichnerowicz conformal approach 

The process of solving of the constraint equations can be significantly simplified 
by the conformal York-Lichnerowicz decomposition [H] , CO : 

7y = '0''%, (6) 
= k,,^~\ (7) 

The new metric 7^ is called the conformal (or background) metric and can 
be chosen arbitrarily. This conformal tranformation preserves the shape of the 
momentum constraints; 

V,{K^^ -rKD^O, (8) 
and the hamiltonian constraint can be rewritten as 

V^V, = - l{K^JK'' - {Klf)i^-'. (9) 

The covariant derivative V and the Ricci scalar R correspond to the background 
metric . We will also impose the maximal slicing condition 

K :^Kl^ 0. (10) 

This equality can be understood as a gauge fixing condition. It fixes the lapse 
function N and is convenient in avoiding coordinate singularities during the 
time evolution of the system. It also simplifies the constraint equations which 
finally take the form 

Vjk^^ = 0, (11) 
= - \k.jk''^p-'. (12) 

2.3 Axially symmetric ansatz 

We would like to solve numerically the constraint equations for the boosted Kerr 
black hole with momentum directed along angular momentum. Therefore it is 
desirable to choose a Weyl metric as the background metric. The Weyl metric 
is axially symmetric and its line element in the spherical coordinates has the 
form [E], [TT] 

%dx'dx^ = exp{~q{r, 9)){dr^ + r^dO^) + sin^ ed(t>^ . (13) 
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Let us recall that the Kerr initial data can be written (in terms of quasi-isotropic 
radial coordinate r) as |3] 

exp(-q(r,e)) = ^, (14) 

T. = R^ + a^cos^e, (15) 

A = i?2-2mi^ + a^ (16) 

E = r + rnH , (17) 

r 



, = (18) 

where m, a are Kerr parameters and the conformal factor ^{t^ 9) is given by 

^\r,e)^^exp{q{r,e)). (19) 

We will fix the exp{—q{r, 9)) factor in (|13|l to be that from the Kerr metric (|14|) . 
The external curvature in the axially symmetric case with maximal slicing can 
be written in the form 

/ll(r,e) /i2(r,9)sine fl3{r,9) 

iK^,)=\ ■■■ h2{r.e) iSp^ I, (20) 

(/ll(r,9) + /22(r,9))sm^e 
exp{—q{r,6)) 

where i,j = r,d,(f> and Kij ~ Kji. The momentum constraint equations Hll|l 
reduce then to 

rdrirfn) + tHhi + f22)drq - ^,((1 - y^)rh2) = 0, (21) 

9.((1 - y'yfu) - ^(/ii + f22)dyq - dy{{l - y')f22) + yhi = 0, (22) 

drhz + def2i = 0. (23) 



The axial symmetry allows one to solve these equations by the method of sep- 
arating variables. The functions /13 and /23 are present in the third equation 
only. This equation can be easily integrated. The result is 

/i3 - -deZ{r,e), (24) 
/23 = drZ{r,e), (25) 

where Z(r, 6) is an arbitrary function. If we put 

hi - /l2 =/22 - 0, (26) 

ry/ n\ / 3 /, o /,^ ma^ sin'' 6* cos(6') 

Z{r,e) = TOa(cos'^ 6* - 3 cos 6*) —, (27) 

then we recover the full Kerr initial data. 
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3 Numerical approach 



3.1 Momentum constraints 

The background metric becomes flat for q{r, 0) = 0. This is the case when 
the angular momentum of the Kerr black hole is equal to zero and we get the 
Schwarzschild solution. Thus for a = we have the fiat background metric and 
then the momentum constraints can be solved analytically (assuming maximal 
slicing and axial symmetry). The solution reads as [SI 

Krr = ^dlW{r,e), (28) 

KrS = -^dydrW{r,e), (29) 
rsmS 

Kee = -^mrdrW{t,e)) + -{ydyW{r,e)-W{r,e))), (30) 
sm 9 r 

Kr^ = ^dyZ{r,e), (31) 
Ks4, = -^drZ{r,e), (32) 

where y = cos 6 and Z are arbitrary functions. The component K^^ can 
be calculated from the maximal slicing condition Kl — 0. These formulae can 
be obtained as a special case of Dain-Friedrich ^ conformally flat initial data 
that have been written in the Newman-Penrose formalism, or more simply by a 
direct integration of the momentum constraint equations (|21I23|) . 

Our idea is to obtain the momentum constraint solutions as a deformation 
of the above flat background formulae with a as the deformation parameter. 
Thus we look for solutions in the form 

k^j (a) = (a = 0) + Sk^j (a) . (33) 

where Kij{a — 0) are given above and 5Kij{a) will be found numerically from 
(|20() . 121I22|I with appropriate boundary conditions which will be discussed be- 
low. 



3.2 Hamiltonian constraint 

Numerical methods are necessary in order to solve the hamiltonian constraint in 
both flat and Weyl background metric. This constraint is a quasilinear elliptic 
equation and can be solved by the standard Newton method. The equation has 
the form 

= F{r,e,i;), (34) 
F{r,e,^) - f{r,e)i,-\ (35) 

where L is a linear operator and F, f are arbitrary functions. Let ipn be the n-th 
approximation in the Newton sequence. Then the following reccurence formula 
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can be easily obtained 



-F^{r,9,1pn-l){i>n - V'n-l)- 

Here Sipn is a correction to the n-th approximation and denotes the partial 
derivative of F with respect to ip- It is important that the conformal factor ip 
p9|l is the analytical solution of ^ for Kerr initial data. Therefore it is an 
excellent example for testing of numerical methods and can serve as the 0-th 
order approximation to more complicated cases. 



3.3 Bowen-York boundary conditions 

Bowen and York demand the metric 7y and the external curvature Kij to be 
invariant (up to a sign) under inversion through a sphere of radius & @]. 
This transformation is given in spherical coordinates by 



r 



(36) 



The metric 7^ is invariant if the background metric jij is flat and if the con- 
formal factor ilj{r,9, (j)) satisfies 



tp{r,9,(j)) = -ip{r, 
r 



(37) 



The same holds true for Weyl background metric if the function g(r, 9) is 
invariant under the inversion transformation. This is the case for the Kerr metric 
with b given in H18|l . The invariance of the extrinsic curvature is equivalent to 
the following transformation rules 



(38) 



We construct a family of invariant solutions of momentum constraints Hll|) for 
the flat background metric by putting the generating function W{r, 9) in for- 
mulae (|28I3UII equal to 



krr{r,9, 


0) = ±^krr{f, 


kre{r,9, 




kee{r,9, 


4>) = ±^kee{f, 


ke^{r,9, 


<^) = T^^ke^if, 


K<f,4,{r, 9, 





■5 U2 

dyW{r, 9) = ~^P{r ± — ) sin^ eg{9), 



(39) 



where P is a parameter and g{9) an arbitrary function. Let us note that for 
g{9) = 1 a Bowen-York solution is usually written in cartesian coordinates in 
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the form 



= — {Pinj+Pjni~{r]ij~ninj)P^nk)T-^{Pinj+Pjni~{r]ij~hnin 

(40) 

where rii is the unit normal to a sphere r = const. In this case Pi is the total 
linear momentum 

P^^^j K^,<fS\ (41) 

directed along z-axis. The appropriate formulae for the Weyl background metric 
result from H39|l plus some corrections (which become zero for a = 0) that 
satisfy symmetry conditions H38|) . We get the following boundary conditions for 
corrections to the extrinsic curvature functions 122|) 

Wii(6,0) + i,5/ii(6,0) = O, 
'5/12(6,0) = 0, 
limr^oo <5/ii(r, 9) = lim^^oo Sfi2{r, 9) = 0. 

The function /22 can be chosen arbitrary. It is reasonable to keep the same 
relation between /n and /22 as in the flat background case, that is 

/ll + ^/22 = 0. 

The boudary conditions for the conformal factor ip result from ^'^7\l (and asymp- 
totical flatness) and read as 

dMb,e) + j-^i;{b,e) = Q, 

limr^oo i'ir, 9) = 1. 
3.4 Puncture boundary conditions 

In the Bowen-York approach each slice consists of two isomorphic ends (the 
inversion transformation (|36|l defines an isometry of the physical metric) sepa- 
rated by a sphere of radius b and we can limit ourselves to one end without any 
singularities. The puncture method takes care of singularities in another way. 
In the flat background case we start with the curvature generated by 

dyW{r, 9) = ~^Prsm^9g{9), (42) 
and continue with the conformal factor of the form 

777 

= 1 + - + <5V, (43) 
Ir 

where Sif) is assumed to be nonsingular on the whole initial slice. Inserting H43|l 
into (O it can be easily seen that the admissible singularities (at r = 0) in 
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curvature functions are of type at most r~^. That is why the term of order 
from had to be excluded in (gU . 

The above formulae can be generalized to the Kerr background case in the 
following manner. The flat curvature functions are completed by (a-dependent) 
corrections in order to satisfy momentum constraint equations (|22|l . These 
corrections are expected to be nonsingular. The form of the conformal factor 
is changed to 

V'pnnc + S'ip, (44) 

where is given in p9|l . The correction 6ip must also be nonsingular. This 
becomes clear if we compare (|43|) with the asymptotic expansions of ip near 
r — and r = oo 

^ir^c^,e) = i + ^ + l^ + o{^). 

Let us note that the boundary conditions for r — s- co are the same as in the 
Bowen-York case, for instance 

lim S^{r,e) = 0. 



4 Numerical results 
4.1 Numerical techniques 

The first step is to find the mapping of the radial coordinate onto some compact 
interval. In the Bowen-York case we have used a new coordinate x defined by 




which mapps the interval [6, oo) onto [0,1]. Similarly, for puncture boundary 
conditions we have put 

r 



In both cases we have to solve elliptic partial differential equations on the rect- 
angle [[0,1], [—1,1]] with mixed boundary conditions. These equations for ex- 
trinsic curvature and conformal factor have been solved on a dense lattice with 
5000x200 points with the help of excellent numerical algorithms for sparse linear 
systems: MUMPS [EI, UMFPACK [H] and HYPRE The precision was 
tested on the conformal factor of the Kerr solution (|19|l . and for instance the 
mass parameter was recovered with the relative error less than 10"'*. 

We have achieved still better precision using finite elements method (FEM) . 
This method is based on triangulations which can be fitted to a shape of suc- 
cesive approximations of the solution. Especially powerfuU has turned out to 
be the program FreeFem-|--|- [S] which is an implementation of the special lan- 
guage dedicated to the finite element method. Wc have used it extensively as 
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it enables one to solve easily PDE problems, both elliptic and time dependent. 
The relative numerical precision has been improved up to 10~^. 

4.2 Detailed results 

Our boundary conditions assume the slices to be asymptotically flat. The ADM 
mass is defined as 

mADM = VE^^Pm, (45) 
where E is the total energy of the system 

E = — -[ (fS'\/i'ilj^\im2r{i:{r,e)~l), (46) 

27r Jg^ r-^oo 

and Pi the total momentum defined in (|41|l . The ADM mass is a parameter 
which shows how quickly a slice becomes flat at infinity. The minimal surface 
is in turn a closed surface which locally has a minimal area (in the riemannian 
geometry) . The following inequality has been proved for slices with the vanishing 
extrinsic curvature dl 



Here m is the ADM mass and S the area of the outermost minimal surface. 
This is the so-called Riemannian version of the Penrose inequality JSIj |19) . 
The minimal surface can be found from the condition 

V,n'' = 0, (48) 

where tt,' is the unit vector normal to the surface. The Penrose inequality must 
be reformulated for slices with non-vanishing extrinsic curvature. The minimal 
surface is then replaced by an apparent horizon by which is understood a closed 
two-dimensional surface which obeys one of the two equations |2()) 

6l± = V,n' T K.jn'n^ ^ 0, (49) 

where ± correspond to the past and future apparent horizons respectively, and 
0± are called optical scalars. In the axially symmetric case these apparent 
horizons satisfy the differential equations |lj 



„3 ^ 

'i;^ _^ qe_ j_ nr^i a\ _ ^2/4Vv - ±1 -I- 1' 



+ ^(^-f +cote)-r2(M._£.^3, 



+roi^ - ^ + cote) - r^{^ - f + ^) = (50) 
with the following boundary conditions for a function r{0) at ~ and = n 

rg\g^O = rg\g=T, ^ 0. (51) 
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The Penrose conjecture for a slice with a non- vanishing extrinsic curvature has 
the same form H47II as before, but S is the area of the outermost apparent horizon 
to the future instead of the minimal surface. The above formulation is one of 
a number of possible wordings of the hypothesis - see |20j for an extensive 
discussion of various formulations of the Penrose conjecture. In fact all the 
versions described in [20] have been checked. Since the results differ very little 
and the hypotheses always stand the tests, we report here only data concerning 
the Penrose inequality described above. This inequality is not proved as yet 
although there are some partial results j > 112] i [2Sj j 123 > jlH] and schemes 
for the proof [20], |23- It is widely believed that if the Penrose inequality is not 
true then the cosmic censorship will be broken. 
Our numerical program was the following one 

• for given values of Kerr parameters m, a solve the momentum constraints, 

• calculate the conformal factor tjj for several values of the black hole mo- 
mentum p = VPiP^ (EH), 

• for each P find the ADM mass, an apparent horizon and its area 



which according to the Penrose inequality should satisfy the inequality 
ex < 1. 

The numerical results are presented on three figures (Fig. . Each of them 
shows how eK depends on momentum in case of Bowen-York and puncture 
boundary conditions. 

Let us note that this dependence is quite different for the Bowen-York and 
the puncture boundary conditions. In the Bowen-York case the coefficient eK 
reaches some maximum, that is evident especially for big values of a. 

There is also another inequality which is expected to hold for every axially 
symmetric, asymptotically flat initial data. Similarly to the Penrose inequality 
we can define a coefficient ea 



where J is the total angular momentum. It is important to note that in our 
case 




graphically show the momentum dependence of the coefficient 





(52) 



J = ma, 



(53) 
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Figure 4: The momentum dependence of eA for a = 0.5 and m = 2. 

even if the total momentum is not zero. The axially symmetric counterpart of 
the Penrose inequality can be stated as 

eA < 1, 

and the equality is reached only for Kerr slices j28| . The momentum dependence 
of the coefficient both for Bowen-York and puncture data is shown on Fig. 

EE. 

Let us stress that our results depend only of two pameters a = — and — . 

^ ^ ^ mm 

Therefore we have limited ourselves to m = 2. 

5 Conclusions 

We have generalized Bowen-York and puncture constructions of black hole ini- 
tial data to the case of a boosted Kerr black hole in the axially symmetric case. 
These initial data smoothly depend on the black hole angular momentum and 
for a — coincide with the conformally flat solutions. Our initial data are very 
precise and can serve as the starting point in constructing initial data for binary 
black holes. It would be also interesting to compute the long-term numerical 
evolution of the single boosted Kerr black hole and look for the quasinormal 
modes present in the radiation. We have calculated numerically the apparent 
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horizons and tested the famous Penrose inequaUty together with its more restric- 
tive version vaHd for axially symmetric initial data only. The numerical results 
are in agreement with theoretical predictions and clearly confirm the fact that 
the equality in the axially symmetric version is possible only for non-boosted 
Kerr black hole. 

Our numerical calculations were mostly based on finite element method 
(FEM) techniques. In our opinion these methods lead to more precise results 
than those made on rectangular grids. The strength of FEM lies in the fact that 
the FEM algorithms have been widely tested in other areas of physics especially 
in elasticity, aerodynamics, electrostatics and in time dependent problems in 
hydrodynamics . 

We believe that our results can be generalized in a straightforward manner 
to the case of binary black hole initial data. The work concerning appropriate 
puncture construction is in progress. 
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supported by the Polish State Committee for Scientific Research (KBN) grant 
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